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(S| ; ABSTRACT 

We measure the power spectrum of the galaxy distribution in the ESO Slice Project 
C*") ■ (ESP) galaxy redshift survey. We develope a technique to describe the survey window 

function analytically, and then deconvolve it from the measured power spectrum using 
, a variant of the Lucy method. We test the whole deconvolution procedure on ESP 

mock catalogues drawn from large N-body simulations, and find that it is reliable for 
recovering the correct amplitude and shape of P{k) at k > 0.065 h Mpc -1 . In general, 
the technique is applicable to any survey composed by a collection of circular fields 
' with arbitrary pattern on the sky, as typical of surveys based on fibre spectrographs. 

, The estimated power spectrum has a well-defined power-law shape k n with n ~ —2.2 

■ for k > 0.2 h Mpc -1 , and a smooth bend to a flatter shape (n ~ —1.6) for smaller 

' fc's. The smallest wavenumber, where a meaningful reconstruction can be performed 

*p\ (k ~ 0.06 h Mpc -1 ), does not allow us to explore the range of scales where other power 

Q ■ spectra seem to show a flattening and hints for a turnover. We also find, by direct 

£j | comparison of the Fourier transforms, that the estimate of the two-point correlation 

c/2 . function £(s) is much less sensitive to the effect of a problematic window function 

as that of the ESP, than the power spectrum. Comparison to other surveys shows 
, an excellent agreement with estimates from blue-selected surveys. In particular, the 

• ■ »h . ESP power spectrum is virtually indistinguishable from that of the Durham-UKST 

' survey over the common range of fc's, an indirect confirmation of the quality of the 

deconvolution technique applied. 
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1 INTRODUCTION 

The quantitative characterisation of the galaxy distribution 
is a major aim in the study of the large-scale structure of 
the Universe. During the last 20 years, several surveys of 
galaxy redshifts have shown that galaxies are grouped in 
clusters and superclusters, drawing structures surrounding 
large voids (see e.g. Guzzo 1999 for a review). The power 
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spectrum of the galaxy distribution provides a concise sta- 
tistical description of the observed clustering that, under 
some assumptions on its relation to the mass distribution, 
represents an important test for different structure forma- 
tion scenarios (e.g. Peacock 1997 and references therein). 
Indeed, under the assumption of Gaussian fluctuations, the 
power spectrum totally describes the statistical properties 
of the matter density field (e.g. Peebles 1980) . 

In recent years, several estimates of the galaxy power 
spectrum have been obtained using galaxy samples selected 
at different wavelenghts: radio (Peacock & Nicholson 1991), 
infrared (Feldman, Kaiser & Peacock 1994, Fisher et al. 
1993, Sutherland et al. 1999) and optical (Park et al. 1994, 
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da Costa et al. 1994, Tadros & Efstathiou 1996, Lin et al. 
1996, Hoyle et al. 1999), to mention the most recent ones. 

The ESO Slice Project redshift survey (ESP, Vettolani 
et al. 1997, 1998) is one of the two deepest wide-angle sur- 
veys currently available, inferior only to the larger Las Cam- 
panas Redshift Survey (LCRS, Shectman et al. 1996). Dur- 
ing the last few years, it has produced a number of statis- 
tical results on the properties of optically-selected galaxies, 
as e.g. the luminosity function (Zucca et al. 1997) or the 
correlation function (Guzzo et al. 2000). The geometry of 
the survey (a thin row of circular fields, resulting in an es- 
sentially 2D slice in space) is such that an estimate of the 
power spectrum represents a true challenge. In this paper 
we present the results of a detailed analysis that overcomes 
these difficulties, producing a reliable measure of the power 
spectrum from the ESP redshift data. The technique devel- 
oped here to cope with the specific geometry of the survey is 
potentially interesting also for application to other surveys 
consisting of separate patches on the sky, as could be the 
case, for example, of preliminary sub-samples of the ongo- 
ing SDSS (Margon 1998) and 2dF (Colless 1998) surveys. 

The outline of the paper is as follows. We shall first 
recall the main features of the ESP survey and the sample 
selection (section |^) , then discuss the power spectrum esti- 
mator adopted for the analysis (section ^) and the numeri- 
cal tests performed in order to check its validity range (sec- 
tion^). We shall then present the estimated power spectrum 
(section 5 ) and its consistency with the correlation function 
(section 6), and then discuss it in comparison to the results 
from other surveys (section Q). Section ^ summarises the 
results obtained, drawing some conclusions. 



2 THE ESO SLICE PROJECT 

The ESO Slice Project galaxy redshift survey (ESP, Vet- 
tolani et al. 1997, 1998) was constructed between 1993 and 
1996 to fill the gap that existed at the time between shal- 
low, wide angle surveys as the CfA2, and very deep, one- 
dimensional pencil beams. The survey was designed in or- 
der to allow the sampling of volumes larger than the max- 
imum sizes of known structures and an unbiased estimate 
of the luminosity function of field galaxies to very faint ab- 
solute magnitudes. The survey and the data catalogue are 
described in detail in Vettolani et al. (1997, 1998). Here we 
limit ourselves to a summary of the main features relevant 
for the present analysis. 

The ESP survey (see Figures [j] and ^|) extends over a strip of 
a x 8 = 22° x 1°, plus a nearby area of 5° x 1°, five degrees 
west of the main strip, in the South Galactic Pole region 
(22 h 30 m < a < 01 h 20 m , at a mean declination of -40°15' 
(1950)). This region was covered with a regular grid of adja- 
cent circular fields, with a diameter of 32 arcmin each, corre- 
sponding to the field of view of the multifibre spectrograph 
OPTOPUS (Avila et al. 1989) at the ESO 3.6 m telescope. 
The total solid angle covered by the survey is 23.2 square de- 
grees and its position on the sky was chosen in order to min- 
imize galactic absorption (—75° ^ b 11 S~ —60°). The target 
objects, with a limiting magnitude bj < 19.4, were selected 
from the Edinburgh-Durham Southern Galaxy Catalogue 
(EDSGC, Heydon-Dumbleton et al. 1989). A total of 4044 
objects were observed, corresponding to ~ 90% of the par- 



ent photometric sample and selected to be a random subset 
of the total catalogue with respect to both magnitude and 
surface brightness. The total number of confirmed galaxies 
with reliable redshift measurement is 3342, while 493 objects 
turned out to be stars and 1 object is a quasar at redshift 
z ~ 1.174. No redshift measurement could be obtained for 
the remaining 208 spectra. As discussed in Vettolani et al. 
(1998), the magnitude distribution of the missed galaxies is 
consistent with a random extraction of the parent popula- 
tion. About half of the ESP galaxies present spectra with 
emission lines. Particular attention was paid to the redshift 
quality and several checks were applied to the data, using 
1) multiple observations of ~ 200 galaxies, 2) ~ 750 galaxies 
for which the redshift from both absorption and emission 
line is available (Vettolani et al. 1998, Cappi et al. 1998). 
More details about the data reduction and sample complete- 
ness are reported in Vettolani et al. (1997, 1998). 

Given the magnitude-limited nature of the survey, the 
computation of a clustering statistics like the power spec- 
trum requires the knowledge of the selection function. This 
is defined as the expected probability to detect a galaxy at 
a redshift z and can be expressed as 



«(*) = 



max[L i 



!(*)] 



(f>(L)dL 



f£°<KL)dL 



(1) 



where cj>(L) is the luminosity function, L\ is the minimum 
luminosity of the sample and L m i n (z) is the minimum lu- 
minosity detectable at redshift z, given the sample limiting 
magnitude. 

In the ESP survey the minimum luminosity corresponds 
to an absolute magnitude Mbj,i = —12.4 + 5 log h (h is the 
Hubble constant in units of 100 kms~ 1 Mpc~ 1 ). L m i n (z) is 
the luminosity of a galaxy at redshift z with an apparent 
magnitude equal to the apparent magnitude limit b.j = 19.4. 
The corresponding absolute magnitude is given by 



M b 



25 + 5 log D L (z) +K(z), 



(2) 



where Dl is the luminosity distance in Mpc and K(z) is 
the K-correction. Dl is given by the Mattig formula (1958), 
which depends on the assumed cosmological model. For all 
ESP computations we assume a flat universe with o = 1 
and A = 0. Before proceeding to the computation of lu- 
minosity distances, we have converted the observed helio- 
centric redshifts in the catalogue to the Cosmic Microwave 
Background (CMB) rest frame using a standard procedure, 
as described in Carretti (1999). The luminosity distance is 
then given by 



D L {z) 



2c 



(1- 



1 



(3) 



The K-correction is a function of redshift and morpho- 
logical type, but the latter is not directly available for ESP 
galaxies. Following Zucca et al. (1997), we use an average K- 
correction, weighted over the expected morphological mix- 
ture at each z. See Zucca et al. (1997, cfr. their figure 1) 
for the details of this computation. A recent principal com- 
ponent analysis of the spectra (Scaramella, priv. comm.) 
confirms the adequacy of this mean correction. 

The luminosity function is such that (f>(L)dL gives the 
density of galaxies with luminosity L G [L, L + dL[. The 
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01 h 20 m 



22 h 30 r 



Figure 1. The area covered by the ESP survey on the sky consists of a set of 107 circular fields of 16' radius. As shown in this figure, 
they are arranged into 2 parallel rows and draw two thin slices over the celestial sphere of about 22° X 1° and 5° X 1° respectively, 
separated by ~ 5° (from Vettolani et al. 1997). 



ESP luminosity function is well approximated by a Schechter 
(1976) function (Zucca et al. 1997) 

*(L)dL = *-(A)V^d(A) (4) 

with best fit parameters a = —1.22, = —19.61 + 5 log ft 
and (j>* — 0.020ft 3 Mpc -3 . In reality, as shown by Zucca et al. 
(1997), for Mb j > —16 + 5 log ft the faint end steepens with 
respect to the Schechter form and the overall shape is better 
described by adding an extra power law. Nevertheless, this 
is relevant only for the very local part of the sample and a 
description of the selection function using a simple Schechter 
fit is perfectly adequate for our purposes. 

Another quantity to be taken into account for clustering 
analyses is the redshift completeness of the 107 fields, as not 
all target galaxies at the photometric limit were succesfully 
measured. This can be expressed as (Vettolani et al. 1998) 

Nz 

Nt-Ns- 0.122N NO ' 
where, for each field, Nt is the total number of objects in the 
photometric catalogue, Nz is the number of reliable galaxy 
redshifts, Nno is the number of not observed objects, N$ is 
the number of stars and 0.122 is the fraction of stars in the 
spectroscopic sample. In Figure ^ we plot the completeness 
values for each field. Field numbers < 100 denote fields in 
the northern row, while the others refer to the southern one. 

The power spectrum analysis has been performed on 
both volume-limited and magnitude-limited subsamples of 
the survey. Volume-limited samples include all galaxies in- 
trinsically more luminous than a given absolute magnitude 
Mii m and within the maximum redshift z max at which such 
magnitude can still be detected within the survey apparent 
magnitude limit. In such a case, the expected mean density 
of galaxies does not vary with distance. Magnitude-limited 
catalogues, by definition, are simply subsets of all galaxies in 
the survey to a given apparent magnitude, possibly with the 
addition of an upper distance cut z max above which the value 
of the selection function becomes too small. Magnitude- 
limited samples contain more objects, but the mean ensem- 
ble properties (as e.g. the galaxy luminosity distribution) 
vary with distance. We extract from the ESP survey two 
magnitude-limited samples with different z max limit, plus 



Table 1. Parameters of the samples extracted from the ESP sur- 
vey: z m ax is the maximum redshift, D max the maximum comoving 
distance in h~ 1 Mpc unit, Mu m the absolute magnitude limit for 
the volume— limited sample (we omit the 5 log h term) and N is 
the galaxy number. 



Sample 



(ft" 1 Mpc) 



Mun 



N 



ESPm523 0.20 523 
ESPm633 0.25 633 
ESP523 0.20 523 



-20.1 



3092 
3306 
481 



one volume-limited sample with Mn m < —20.1 + 5 log ft. 
(For simplicity, we shall omit hereafter the 5 log ft term). 

For the estimate of the power spectrum, comoving dis- 
tances are computed for each galaxy as D c (z) — Dl(z)/(1 + 
z). The uncertainty introduced in D c because of our igno- 
rance of the correct cosmological model amounts to less than 
5% for a typical redshift z = 0.20, when the value of Q is 
changed from 1 to 0.2. 



3 POWER SPECTRUM ESTIMATOR 

The galaxy power spectrum can be defined as 



P(k) 



(,(x)e 1 x dx, 



(6) 



where £(x) is the two-point correlation function, x and k are 
the comoving position and wavenumber vectors respectively, 
while x = |x| and k — |k|. Under the hypothesis of homo- 
geneity and isotropy, P(k) and £(x) are functions only of k 
and x respectively. By definition the two-point correlation 
function can be also written 



P{k) oc \5{k)\ 2 



(7) 



where <5(k) is the Fourier transform of the density contrast 
of the galaxies. 

In this paper we follow the Fourier notation 
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-40.8° < 6 < -39.7° 



Figure 2. The galaxy distribution in the ESP redshift survey. 



/>)= j /(x)e-*- x dx, (8) 
/(x) = ^y./(k)e lkx dk. (9) 

To compute the power spectrum of galaxy density fluc- 
tuations from the observed galaxy distribution, we use a tra- 
ditional Fourier method (cfr. Carretti 1999 for details), as 
developed by several authors (e. g. cfr. Peebles 1980, Fisher 
et al. 1993, Feldman et al. 1994, Park et al. 1994, Lin et 
al. 1996). We also apply a correction (Tegmark et al. 1998), 



that accounts for our ignorance on the true value of the mean 
density of galaxies (Peacock & Nicholson 1991). 

Given a sample of N galaxies of positions Xj and weights 
Wj , an estimate of the Fourier transform of density contrast 
is given by 

<5(k) = -J- w ^ 1 - WW, (10) 

where V is the volume of the sample and VK(k) is the Fourier 
transform of the survey window function (hereafter a ~ will 
denote the quantities estimated from the data) . The window 
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Figure 3. The redshift completeness within the ESP fields, i.e. the fraction of galaxies with measured redshift with respect to the total 
number of galaxies in the photometric sample in each field. Field numbers are as reported in the catalogue (Vettolani et al., 1998). The 
two panels are for the fields in the northern (a) and southern (b) rows respectively. 



function W(x) is 1 within the volume covered by the sample 
and elsewhere, so it can be described as an ensemble of 
107 cones. This geometry allows us to obtain analytically the 
Fourier transform VF(k) as the sum of the Fourier transform 
of each cone. In equation [uj each galaxy contributes with 
some weight Wj. In a volume-limited catalogue all galax- 
ies have equal weight, i.e. Wj = 1. In a magnitude-limited 
catalogue the expected galaxy density decreases with the 
distance according to the selection function. Thus, the sim- 
plest form for the weight for a galaxy is given by the inverse 
of the selection function at its redshift Zj 

1 



(11) 



If the catalogue completeness is C < 1, the previous weight 
should be modified as 

for volume-limited catalogues, and as 
1 



(12) 



(13) 



3 s{z d )c{^y 

for magnitude-limited catalogues. C(xj) is the completeness 
of the sample at the position of the j th galaxy. 

Our adopted power spectrum estimator is defined with 

respect to <5(k) by the following equation (Tegmark et al. 
1998) 



A(k) = 



where 



l^(k)| 2 -6(k) 
A(k) ' 



(14) 



A(k) = 1 - 



W(k) 



W(0) 



V 



(15) 



accounts for our ignorance of the mean galaxy density, while 



6(k) = 



V 



3=1 



W(k) 



W(0) 



(16) 



is the shot noise correction due to the finite size of the sam- 
ple. 

The observed power spectrum estimated by equation |l4| 
is related to the true power spectrum P(k) by 



^c(k) = 

where 



(2 7 r) 3 J 4(k) 



P(fe')0(k,k')dk', 



W(0) 



(17) 



(18) 



For wavenumbers k such that |W"(k)| <<iC 11^(0)1 this equa- 
tion reduces to the convolution between P(k) and ^(k)] 2 . 

To describe the convolved power spectrum we choose to 
average -P c (k) over all directions 



P c {k) = (P c (k)> 



4tt 



Pc(k)<m k 



k' P(k')x(k,k')dk', 



(19) 
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where fife is the sphere defined by wavenumbers of amplitude 
k. The kernel of this integral equation is given by 



X (k,k') 



where 



2(2tt) 4 V 



V(k,k')dn A rfn fe /, 



V>(k,k') = 



I W(k - k') - W(-k') W(k)/W(0) | 
1 - I W(k)/W(0)\ 2 



(20) 



(21) 



and f2 fc / is the sphere defined by wavenumbers of amplitude 
k'. 

The Fourier transform of the window function has been 
analytically computed as the sum of the Fourier transforms 
of all the 107 cones. In reality, the cones are slightly over- 
lapped but the small common volume (2.85%) allows us to 
make the assumption of disjoined cones. 

The small width of one ESP cone allows us to ana- 
lytically compute its Fourier transform. Let r be the cone 
height and A9 < 1 rad its width {AO = 16' = 0.00465 rad). 
In the simple case of a cone centered on the z axis, the 
Fourier transform is 



w (k) 



dr r 

JQ 



d9 sinr?e 



(22) 



Taking into account the small value of A9, the integrand 
can be approximated to first order in 8, resulting in 



W c (k) = 



dr r 
o jo 



d9 sin 9 e 



2tt[1 - cos(A0)] / drr 2 e- lhrcoBa . (23) 



This expression depends on r , A8 and cos a, where a is the 
angle between k and the z axis. For a generic cone along the 
direction (9 , <f> ) a rotation can be applied in order to bring 
the cone on the z axis. So, the Fourier transform is 



w c (k) 



dr r 
o J si 



dfle~ 



(8o,*o) 



» 2tt[1 - cos(A0)] / drr 2 e- lkrcoB ' 1 , (24) 
Jo 

where tt($ ,$ ) is the solid angle of the cone centered on 
(9 , <j)o), and 7 is the angle between the wave number k and 
the direction (9 ,(j> ). By solving the integral one gets 



w c (k) 



2 ^[l-cos ( A0)] x 



fccos7 



2 -i/Er n cos 7 .2t"o6 

r n e — z 



-ikr D cos 7 



-ikr cos "v 



k cos 7 



(fccos 7) 2 



(25) 



Finally, the Fourier transform of the whole ESP window 
function is given by the sum of the Fourier transform of the 
cones normalized to the true volume of the survey to account 
for the overlapping zone 



W(k) = 



V 



107 

v: ^ 

i=l 



(26) 



where i is the cone index, V c is the volume of one cone and 
V is the true survey volume 



V = 107 14(1-/3), 



(27) 



10 




0.1 

k (h Mptr 1 ) 



Figure 4. The power spectrum of the ESP window function (here 
we use this term for the real-space survey selection window, but 
in other papers the same is used for its Fourier transform di- 
rectly). The filled line is the spherically-averaged function, com- 
puted analytically using the machinery described in the text. The 
filled points give the same quantity, but computed numerically 
through a simple Montecarlo simulation. The figure shows also 
the three components of |W(fc)| 2 in k-space. Note how broad is 
the function along the direction 2, which has been chosen to be 
essentially perpendicular to the ESP main plane, evidencing its 
extreme anisotropy. 



which accounts for the total volume fraction (3 = 0.0285 lost 
in the cone overlaps. 

To check the reliability of our analytic approximation, 
we perform a numerical Fourier computation. We sample 
the survey volume on a regular grid by assigning 1 to the 
grid points inside the window function and outside. We 
then perform an FFT. This numerical computation is lim- 
ited by the finite size of the grid cells, but avoids the overlap- 
ping zones and considers the true window function. Figure ^ 
(filled points and solid line) compares both the analytical 
and numerical estimates of the window function power spec- 
trum averaged over spherical shells. The difference is less 
than 5%. 

The strongly anisotropic geometry of the ESP survey 
(see Figure ^J) introduces important convolution effects be- 
tween the survey window function and the galaxy distribu- 
tion. To clean the observed power spectrum for these effects, 
we have adopted Lucy's deconvolution method (Lucy 1974; 
see also Baugh & Efstathiou 1993 and Lin et al. 1996, for 
a discussion about its application to power spectrum esti- 
mates) . 

The Lucy technique is a general method to estimate 
the frequency distribution V( r ?) °f a quantity rj, when we 
know the frequency distribution <j>(y) of a second quantity 
y, related to r\ by 



4>{y) 



ip(r])U(y\T]) drj, 



(28) 
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where H(y\r)) dy is the probability that y' G [y, y+dy[ when 
rj = rj. The probability H(y\n) must be known and the fre- 
quency distribution 4>{y) is the observed one. 

The solution of equation ^ can be obtained by an iter- 
ative procedure. Let Q(ri\y) dr\ be the probability that rj £ 
[rj, ri+dri[ when y — y. The probability that y G [y, y + dy[ 
and rj G [rj, n + drj[ can be written as (f>(y)dy Q(n\y)dn 
and ip(n)dnll(y\n)dy. From these two expressions and equa- 
tion we obtain 

i>{rj)TL{y\r)) 



® VV J ij)(r))Il(y\ri) dn' 
which provides the identity 

i>(v) = / 4>{y)Q{v\y)dy. 



(29) 



(30) 



The latter equation cannot be solved directly, since 
Q(?l\y) depends on the unknown ip(n) as well. Given a fidu- 
cial model for i()(n) and the known probability II(y|77), equa- 
tion by provides an estimate for Q(n\y). This and the iden- 
tity pOj allows us to compute an improved estimate for ip(n). 
The process can then be repeated until convergency. In our 
specific case the equation to be solved is eq. [19l where 
k' 2 x(k,k') plays the role of the probability TL{y\rf). If we 
sample k on logarithmic intervals the convolution integral 
becomes 



Pc(fc) = / fc' 3 P(fc')x(M')MlO) d(log 10 fc') 



(31) 



and an iterative scheme for the deconvolved spectrum can 
be written as 



=P m (ki) 



£. [Pc(k j )/P c m (k j )] x (k h k i 



T,jX(kj,ki) 
where 



(32) 



(33) 



and A = (log 10 k i+ i — log 10 ki) is the logarithmic interval, 
while P m denotes the m th estimate of the spectrum. 

One problem with the Lucy method is that of produc- 
ing a noiser and noiser solution as the iteration converges. 
To avoid this, Lucy suggests to stop the iteration after the 
first few steps. This is quite arbitrary and we prefer to fol- 
low Baugh & Efstathiou (1993; see also Lin et al. 1996) in 
applying a smoothing procedure at each step 

P m {h) =0.25P m (fc 1 :_ 1 )+0.5P m (fc ! ) + 0.25P m (fc 1+1 ). (34) 

We use P°(ki) = constant as initial guess for the power 
spectrum, but we checked that the solution is independent 
of the shape of P°(ki). One consequence of this smoothing 
is that some degree of correlation is introduced among the 
bins of P(k). 

The importance of the convolution effects on different 
scales can be estimated by plotting the integrand of eq. ( pjlj ) 
as a function of k', for different values of k (Figure |E|). If 
the window was a large and regular sample of the Universe, 
the plots would be sharply peaked at k — k' , as it actually 
happens for large values of k (small scales). On the other 
hand, for small values of k, i.e. for spatial wavelengths com- 
parable to the typical scales of the window (which are quite 
small, due to the strongly anisotropic shape), the true power 
is spread over a wide range of wavenumbers. 
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Figure 5. Behaviour of the integrand of the convolu- 
tion equation normalized with respect to the convolved 
power spectrum P(k) for some k values. We plot the 
quantity f(k') = k' 3 P(k') X {k, k 1 ) ln(10)A/P(fc) for (k = 
0.032, 0.057, 0.1, 0.18, 0.32/i Mpc" 1 ), considering the kernel rel- 
ative to the geometry of the 523 /i — 1 Mpc sample. The power 
spectrum P(k') of a CDM model with Q = 0.4 and V = 0.2 has 
been used. 



4 NUMERICAL TESTS 

We test the whole procedure for estimating the power spec- 
trum through iV-body simulations that we have run assum- 
ing some cosmological models (Carretti 1999). The results 
of these simulations can be considered as a Universe, from 
which we can extract mock catalogues with the same fea- 
tures of the ESP survey (geometry, galaxy density, field 
completeness, selection function). We then apply to such 
mock catalogues the whole power spectrum estimate pro- 
cedure (convolved power spectrum estimator and deconvo- 
lution tecnique) and we compare the result with the true 
power spectrum obtained from the whole set of particles. 

The simulations were performed on a Cray T3E 
at CINECA supercompunting center (Bologna) using a 
Particle-Mesh (PM) code (Carretti & Messina 1999) and 
adopting two cosmological models: an unbiased Standard 
Cold Dark Matter (SCDM: Q c = 1, h = 0.5, a 8 = 1) and 
an unbiased Open Cold Dark Matter with shape parameter 
T = 0.2 (OCDM: D. = 0.4, h = 0.5, a 8 = 1). They were 
run with a box size of 700ft" 1 Mpc, 512 3 grid points and 
512 3 particles, in order to reproduce a volume which can 
contain all catalogues selected for the analysis (max depth 
633/i -1 Mpc) and to select a realistic number of mock galax- 
ies for the magnitude-limited catalogues. From each simula- 
tion box, we randomly choose a particle as origin and extract 
sets of particles with the same features of the three ESP cat- 
alogues. The magnitude-limited selection is then reproduced 
by simply assigning a weight corresponding to the observed 
selection function. 

From each simulation and for each ESP subsample we 
construct 50 independent mock catalogues. The average 
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Figure 6. Deconvolved power spectra from 50 mock ESP catalogues extracted from SCDM (left panels) and OCDM (right panels) 
simulations. The filled squares and the error bars give the ensemble average and standard deviation of the 50 mock samples. The solid 
line is the corresponding power spectrum computed from the whole simulation box using all particles. N denotes the average number 
of particles among the mock catalogues. The three pairs of panels, from top to bottom, refer to the three different kinds of subsamples 
ESPm523, ESP523, and ESPm633, as in the case of the real data. In panel a) we plot also the power spectrum estimate before applying 
the deconvolution procedure (open circles). 
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number of particles over the 50 realisations is set to the 
number of galaxies observed in the corresponding true ESP 
sample. The power spectrum estimator is then applied to 
each of the 50 mock galaxy catalogues, producing an inde- 
pendent estimate of P(k) for that specific model and sample 
geometry. From each set of realisations, a mean P(k) and its 
standard deviation can finally be computed and compared 
to the true power spectrum obtained from all the particles. 
A general result from this exercise is that for k > 

0. 065/i Mpc -1 the systematic power suppression by the win- 
dow function convolution is properly corrected for by our 
procedure, i.e. we are able to fully recover the input P(k). 
In figure H we show the reconstructed power spectra from 
the three sets of mock catalogues, compared to the "true" 
ones, both for the SCDM and OCDM simulations. In panel 
a), in particular, we also show (open circles) the raw power 
spectrum before applying the Lucy deconvolution, to em- 
phasize the dramatic effect of the ESP window function on 
all scales. It is evident that for k > 0.065/i Mpc - the mean 
deconvolved power spectra are a very good reconstruction 
of the original ones. In particular, in the SCDM case where 
the spectrum turnover scale is well sampled by the simula- 
tion box, the technique is able to nicely follow the change 
of shape at small fc's. This is important, because guaran- 
tees that the deconvolution method has enough resolution 
as to follow possible features in the data power spectrum, 
while at the same time recovering the correct amplitude. At 
smaller fe's the error bars explode, and the results become 
meaningless. In the case of the deepest sample, ESPm633, 
the reconstruction shows a small systematic overestimate of 
the amplitude for the SCDM spectrum on the largest scales, 

1. e. the reconstruction algorithm seems to have difficulty in 
following the curvature of the spectrum accurately. This is 
probably due to the rather small value of the selection func- 
tion in the most distant part of the sample, which puts too 
large a weight on the distant objects. Rather than indicating 
a difficulty in the technique, this is probably telling us that 
it is safer to truncate the data at smaller distances, as for 
sample ESPm523. 

Comparing the results from the SCDM and OCDM 
mock catalogues, we have checked that the fractional errors 
for the two cases are quite similar. Not knowing a priori 
the correct cosmological model, rather than choosing one of 
the two models as representative, we prefer to average the 
fractional errors measured from the two models. 

Using the mock catalogues, we can also evaluate directly 
the possible effects of the field incompleteness on the power 
spectrum estimate. Figure ^ compares the mean power spec- 
tra obtained from one set of 50 SCDM mock samples both in 
the ideal case (all fields complete) and when the ESP field- 
to-field incompleteness is introduced and corrected for. It 
is clear that the incompleteness is correctly taken into ac- 
count by the weighting scheme. Error bars (not reported for 
clarity) are also very similar. 
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Figure 7. Test for the effect of redshift incompleteness. The 
power spectrum has been computed for 50 mock samples ex- 
tracted from the SCDM simulation, both in the ideal case of a full 
redshift coverage of a sample as ESPm523, and in the real situa- 
tion, i.e. including the field-to-field incompleteness and correcting 
for it through the weighting scheme. 



spectrum, so that we can now apply it to the three galaxy 
subsamples ESPm523, ESP523 and ESPm633. 

The final results of the computation are shown in fig- 
ure |^ and Table [| The error bars are partially reported only 
for ESPm523 to avoid confusion (the errors are similar for 
the three samples). The three estimates of the power spec- 
trum are well consistent with each other. Given the large 
amplitude of the errors (~ 30% of P(k) for k > 0.15/i 
Mpc" 1 , 50% for k ~ 0.1 and 75% for k ~ 0.065) the small 
differences in the slopes are not significant. In general, we 
can safely say that the power spectrum of ESP galaxies fol- 
lows a power law P(k) oc k n with n ~ —2.2 for k > 0.2h 
Mpc -1 , and n ~ —1.6 for k < 0.2h Mpc -1 . In the range 
0.065 < k < 0.6 h Mpc -1 there is no meaningful difference 
between the three estimates, which are therefore indepen- 
dent of the catalogue type (magnitude- or volume-limited) 
and of the catalogue depth. 

In figure |^ we have also plot, for comparison, redshift- 
and real-space power spectra computed from the simula- 
tions described in section @ (SCDM and F = 0.2 OCDM). 
Note how the former compensate for non linear evolution at 
large k's in this way steepening the slope of P(k) over the 
whole observed range, which makes the global slope closer 
to the observed one. Despite this comparison to models is 
deliberately limited, one can safely say that the data points 
(especially below k = 0.1-0.2 h Mpc -1 ) are in better agree- 
ment with the power spectrum of the V = 0.2 OCDM model. 
This model would reproduce this observation without bias- 
ing (the normalisation adopted by the simulation is as = 1). 



The numerical tests performed have given us an estimate of 
the reliability of our method to reconstruct the true power 
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Figure 8. The final deconvolved estimates of the ESP power spectrum from the three samples. The plot shows also the power spectra 
computed from the two simulations described in the text both in redshift— and in the real— space. Note how redshift distorsion effects 
modify the power spectra increasing the apparent power on large scales and reducing it on small ones. 



Table 2. Results of the power spectrum estimates from the three subsamples of the ESP survey. Errors are estimated from 50 mock 
realisations of the samples, as detailed in the text. 



k (TlMpc 1 ) P{k)ESPm523 lo-ESPm523 P(k)ESPm633 l°",ESPm633 P(k)ESP523 1°"ESP523 



0.065 


20284.1 


15313.7 
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0.075 
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14907.9 
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9786.9 
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0.100 
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1964.7 


638.8 


1983.4 


674.5 
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1187.8 


333.6 


1557.9 


505.0 


1455.4 


500.3 


0.365 


853.9 
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1216.7 
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1042.6 


360.9 


0.422 
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295.7 
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449.0 
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71.5 


151.4 
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Figure 9. Fits with a simple phenomenological form of the convolved and de— convolved P(k) from the ESP523 sample 
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Figure 10. The Fourier transform of the convolved (dashed line) and de-convolved (solid line) estimates of the power spectrum compared 
to the two— point correlation function of the ESP survey (filled circles), estimated for essentially the same volume-limited subsample 
(Guzzo et al. 2000). The Fourier transform of the de— convolved estimate is in very good agreement with the direct measure of £(s). 
This result shows also how the two-point correlation function — for which no kind of correction has been applied in addition to those of 
standard estimators — is substantially insensitive to the effect of the window function. 
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6 CONSISTENCY BETWEEN REAL AND 
FOURIER SPACE 

It is interesting to compare the Fourier transform of the 
ESP power spectrum estimated in this work with the two- 
point correlation function measured independently from the 
same sample (Guzzo et al. 2000). This exercise is a further 
check of the robustness and self-consistency of the estimate 
of P(k). In addition, it is of specific interest to verify the 
effect of the survey geometry /window function in real and 
Fourier space. To simplify the procedure, we have first fitted 
the observed P(k) with a simple analytical form with two 
power laws connected by a smooth turnover (e.g. Peacock 
1997) 

(k/k a ) a 



P(k) 



(35) 



l + (k/k c ) a ~ n 

where k is a normalisation factor, k c gives essentially the 
turnover scale, n is the large-scale primordial index (here 
fixed to n — 1), and a gives the slope for k 3> k c . We have 
used this function to reproduce the global shape of both 
the convolved and de-convolved estimates of P(k) from the 
ESP523 sample. In terms of selection function, this sample 
is the closest to one of the volume-limited samples used by 
Guzzo et al. (2000) to estimate £(s) from the same data. 
Figure shows how this form provides a good description of 
the ESP power spectrum, with the deconvolved one char- 
acterised by k D = 0.080 h Mpc" 1 , k c = 0.062 ft Mpc" 1 , 
a — —2.2 (note that while the slope a is a stable value, 
the turnover scale k c is very poorly constrained, given the 
limited ran ge c overed by the data). 

Figure hoi shows the Fourier transform of the two fits, 
compared to the direct estimate of £(s) by Guzzo et al. 
(2000). Two main comments should be made here. First, our 
"best" estimate of P(k), deconvolved for the ESP window 
function according to our recipe, reproduces rather well the 
observed two-point correlation function (solid line). Note 
how the Fourier transform of the simple direct estimate suf- 
fers from a systematic lack of power as a function of scale 
(dashed line), as we expected from our results on the mock 
samples. The second, more general comment concerns the 
stability of the two-point correlation function. One might 
naively think that the narrowness of the explored volume, 
which gives rise to the window function in Fourier space, 
should affect in a similar way also the estimate of clustering 
by the two-point correlation function. Figure |l(j shows that 
this is not the case. In fact, the points showed here have not 
been subject to any kind of correction (Guzzo et al. 2000) , a 
part from those which are standard in the estimation tech- 
nique to take into account the survey boundaries. Still, they 
seem to sample clustering to the largest available scales in 
a reasonably unbiased way, without basically being affected 
by the survey geometry. 



7 COMPARISON TO OTHER REDSHIFT 
SURVEYS 

In the six panels of figure [ll| we compare the power spec- 
trum for the ESPm523 and ESPm633 samples to a variety 
of results from previous surveys, both selected in the opti- 
cal and infrared (IRAS) bands. In general, there is a good 
level of unanimity among the different surveys concerning 



the slope of P(k) over the range sampled by the ESP es- 
timate. Optically-selected surveys show a good agreement 
also in amplitude, with a possible minor differential biasing 
effect in the case of CfA2-SSRS2-130 (panel a, da Costa 
et al. 1994), which is a volume-limited sample containing 
galaxies brighter than ~ M* — 1.5. The effect of different 
biasing values is more evident in the comparison to IRAS- 
based surveys (IRAS 1.2 Jy, Fisher et al. 1993; QDOT, Feld- 
man et al. 1994; PSCz, Sutherland et al. 1999) in panels e 
and f. 

Particularly relevant is the comparison to the results of 
the Durham/UKST galaxy redshift survey (DUKST, Hoyle 
et al., 1999) (panel d). This survey is selected from the same 
parent photometric catalogue as the ESP (the EDSGC) and 
contains a comparable number of redshifts. However, it is 
less deep (bj 17), while covering a much larger solid angle by 
measuring redshift in a sparse-sampling fashion, picking one 
galaxy in three. This produces a window function which is 
essentially complementary to that of the ESP survey, with a 
good sampling of long wavelengths and a poor description of 
small-scale clustering, which on the contrary is well sampled 
by the ESP 1-in-l redshift measurements. The agreement 
between these two data sets is impressive. This is a further 
confirmation of the quality of the deconvolution procedure 
we have applied to the ESP data, given the rather 3D shape 
of the DUKST volume which makes the window function 
practically negligible for this survey. Significantly more noisy 
is the estimate from the similarly 6j-selected Stromlo-APM 
redshift survey (Tadros & Efstathiou 1996), most probably 
because of the very sparse sampling of this survey and the 
smaller number of galaxies. 

Finally, panel b) shows a comparison with the data from 
the r-selected LCRS (Lin et al. 1996). The power spec- 
trum from this survey has a flatter slope with respect to 
our estimate from the ESP. More in general, it is flatter 
than practically all other power spectra shown in the figure. 
This is somewhat suspicious, as the two-point correlation 
functions agree rather well for ESP, LCRS, Stromlo-APM 
and DUKST (Guzzo 1999), and might be an indication that 
the effect of the window function has not been fully removed 
from the estimated spectrum. 



8 SUMMARY AND CONCLUSIONS 

The main results obtained in this work can be summarised 
as follows. 

• We have developed a technique to properly describe the 
ESP window function analytically, and then deconvolve it 
from the measured power spectrum, to obtain an estimate 
of the galaxy power spectrum. The tests performed on a 
number of mock catalogues drawn from large A r -body simu- 
lations show that the technique is able to recover the correct 
shape of P(k) down to wavenumbers k ~ 0.065 h Mpc -1 . In 
general, this technique for describing the window function 
analytically can be applied to any redshift survey composed 
by circular patches on the sky (e.g. the ongoing 2dF sur- 
vey). In addition to its mathematical elegance, it has some 
computational advantages over the traditional method for 
recovering the survey window function, normally based on 
the generation of large Montecarlo poissonian realisations. 
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Figure 11. Comparison of the ESP P(k) with results from other surveys, as indicated by the labels (see text for references). Error bars 
for the ESP points are reported only partially for clarity. 
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• The final estimates of the ESP power spectrum, ex- 
tracted from three subsamples of the survey, are in good 
agreement within the error bars. The bright volume-limited 
sample does not show a clear difference in amplitude with 
respect to the apparent-magnitude limited ones. This agrees 
with the similar behaviour found for the two-point correla- 
tion function, i.e. a negligible evidence for luminosity seg- 
regation even for limiting absolute magnitudes M&j ~ —20 
(Guzzo et al. 2000). This is only apparently in contrast with 
the results of Park et al. (1994), who found evidence for lu- 
minosity segregation studying the amplitude of the power 
spectrum in the CfA2 survey. In fact, that analysis con- 
centrates on a range of luminosities about 1.5 magnitude 
brighter than M*, which for the CfA2 survey has a value of 
-18.8 (Marzke et al. 1994), i.e. nearly one magnitude fainter 
than for the ESP. This also agrees with the results of Benoist 
et al. (1996), who studied the correlation function for the 
SSRS2 sample, finding negligible signs of luminosity segre- 
gation for M > M* . 

• All three estimates of P(k) show a similar shape, with 
a well defined power-law k n with n ~ —2.2 for k > 0.2 h 
Mpc -1 , and a smooth bend to a flatter shape (n ~ —1.6) 
for smaller fc's. The smallest wavenumber where a meaning- 
ful reconstruction can be performed (k ~ 0.065 h Mpc ), 
does not allow us to explore the range of scales where other 
power spectra seem to show a flattening and hints for a 
turnover. In the framework of CDM models, however, the 
well-sampled steep slope between 0.08 and 0.3 h Mpc -1 
favours a low-T model (P = 0.2), consistently with the most 
recent CMB observation of BOOMERANG/MAXIMA ex- 
periments (Jaffe et al. 2000). 

• We have verified that the two-point correlation function 
£(s) is much less sensitive to the effect of a difficult window 
function as that of the ESP, than the power spectrum. In 
fact, the measured correlation function (without any correc- 
tion), agrees with the Fourier transform of the power spec- 
trum, only after this has been cleaned of the combination by 
the window function. This is an instructive example of how 
these two quantities, despite being mathematically equiva- 
lent, can be significantly different in their practical estimates 
and be very differently affected by the peculiarities of data 
samples. 

• When compared to previous estimates from other sur- 
veys, the ESP power spectrum is virtually indistinguishable 
from that of the Durham-UKST survey over the common 
range of wavenumbers. In particular, between 0.1 and 1 
h Mpc -1 our power spectrum has significantly smaller er- 
ror bars with respect to the DUKST, by virtue of its su- 
perior small-scale sampling. The absence of any systematic 
amplitude difference between these two surveys - both se- 
lected from the EDSGC catalogue, but with complementary 
volume and sampling choices - is an important indirect in- 
dication of the quality of the deconvolution procedure ap- 
plied here, and also of the accuracy of the two independent 
estimates. In this respect, a combination of the Durham- 
UKST and ESP surveys possibly provides the current best 
measure of P(k) for blue-selected galaxies in the full range 
~ 0.03 — 1 h Mpc -1 . It will be very interesting to compare 
these combined results to the power spectrum of the forth- 
coming 2dF redshift survey, which is also selected in the 
same b.j band to virtually the same limiting magnitude than 
the ESP. 
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